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ABSTRACT 

Following up our recent analysis devoted to the impact of non steady accretion on the location of young low-mass stars or 
brown dwarfs in the Herzsprung-Russell diagram, we perform a detailed analysis devoted to the effect of burst accretion on 
the internal structure of low-mass and solar type stars. We find that episodic accretion can produce objects with significantly 
higher central temperatures than the ones of the non accreting counterparts of same mass and age. As a consequence, 
lithium depletion can be severely enhanced in these objects. This provides a natural explanation for the unexpected level 
of lithium depletion observed in young objects for the inferred age of their parent cluster. These results confirm the limited 
reliability of lithium abundance as a criterion for assessing or rejecting cluster membership. They also show that lithium is 
not a reliable age indicator, because its fate strongly depends on the past accretion history of the star. Under the assumption 
that giant planets primarily form in massive disks prone to gravitational instability and thus to accretion burst episodes, the 
same analysis also explains the higher Li depletion observed in planet hosting stars. At last, we show that, depending on 
the burst rate and intensity, accretion outbursts can produce solar mass stars with lower convective envelope masses, at ages 
less than a few tens of Myr, than predicted by standard (non or slowly accreting) pre-main sequence models. This result has 
interesting, although speculative, implications for the recently discovered depletion of refractory elements in the Sun. 

Key words. Stars: formation — Stars: low-mass — Stars: abundances — Accretion, accretion disks 
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1. Introduction 



There is a growing consensus in the star formation commu- 
nity that non steady (episodic) accretion plays a dominant 
role during the formation of low-mass stars (see e.g Enoch 
et al. 2009; Vorobyov 2009; Zhu et al. 2010a and refer- 
ences therein). In a recent paper (Baraffe, Chabrier, Gallardo 
2009, hereafter BCG09), we have suggested that episodic 
accretion provides a viable explanation for the observed lu- 
minosity spread in young cluster Herzsprung-Russell dia- 
grams (HRD). The present follow-up analysis explores in 
more details the effects of episodic accretion on the inter- 
nal structure of young low mass stars (< 1 M Q ), commonly 
used to derive ages of star forming regions and young clus- 
ters. We show that, depending on the accretion history, the 
internal structure of these objects can be strongly affected 
for up to a few tens of Myr (§|2]i. Lithium depletion and, 
for the partly convective stars, the size of the convective en- 
velope, in particular can strongly differ from the standard 
(non accreting) pre-main sequence model predictions (§0. 
In section|U we examine the impact of episodic accretion on 
the observational signatures and show that taking this pro- 
cess into account in the young low-mass object evolution 
provides a consistent explanation for the puzzling observa- 
tions of strong lithium depletion in several low-mass stars 
(LMS) belonging to young clusters (e.g Kenyon et al. 2005; 
Sacco et al. 2007) and in planet-hosting stars (Israelian et al. 
2009), as well as to the recently determined peculiar abun- 
dances of refractory elements in the Sun (Melendez et al. 
2009; Ramirez et al. 2009). 



2. Effect of mass accretion on the internal 
structure 

2.1. Evolutionary models with accretion 

We adopt the same input physics and the same treatment 
of accretion as outlined in BCG09. In standard stellar evo- 
lution calculations, energy conservation equation for a non 
accreting object reads as 



dm), \ dt I 



+ &i 



(1) 



where m is the mass enclosed in a sphere of radius r within 
the object, S the specific entropy and e nuc the local nuclear 
energy generation rate. For very young objects, only deu- 
terium fusion provides a contribution to e nuc [l For accretion 
proceeding at a rate M, time derivatives at fixed mass shell 
must account for the variation of mass with time and of en- 
tropy with accreted mass, i.e. (Sugimoto & Nomoto 1975): 



dS\ (dS\ (dS\ 



(2) 



with q = m/M„. The mass (M-At) accreted during a timestep 
At is assumed to be redistributed over the entire structure, 
with the new mass in a shell of fixed q given by 



m(t + At) = q x (M*(f) + M At). 



(3) 



1 The nuclear energy generation produced by lithium fusion is 
completely negligible. 



2 



Baraffe & Chabrier: Effect of episodic accretion on internal structure and lithium depletion 



In the present calculations, we assume instantaneous 
and uniform redistribution of the extra source of internal en- 
ergy brought by the accreted material. In reality, mass and 
heat redistributions inside the accreting object depend on 
the thermal properties of the accreting material. Proto low- 
mass stars below about 2 M & are expected to be entirely 
convective (Stahler & Palla 2005). For completely convec- 
tive objects, given the short typical convective timescales, 
our assumption of uniform mass and heat redistribution 
should be valid, providing the entropy of accreted matter 
is comparable to or less than the accreting object's inter- 
nal one, so that the infalling material can rapidly thermalize 
with its surroundings (see discussion in Siess & Forestini 
1996; Hartmann et al. 1997). If the aforementioned condi- 
tion is not fulfilled, the accreting matter may not be able 
to penetrate very deep inside the object, possibly leading to 
Rayleigh-Taylor like instabilities in stably radiative regions. 
Siess et al. ( 1 997) explored the effect of mass and heat redis- 
tribution by using a penetration function. Such an approach, 
however, remains highly phenomenological and thus of lim- 
ited reliability. Given the complexity of the problem and 
the present exploratory nature of the effect of episodic ac- 
cretion on young low-mass objects, not mentioning uncer- 
tainties in the accretion processes, we elected to stick to 
the simplest and probably most reasonable in most of the 
presently explored situations treatment, and assumed homo- 
geneous, rapid heat and matter redistribution within the ac- 
creting body. 

Omitting the fraction e( 1 - a) G> £ of the accretion 
shock energy radiated away (see Eqs. (1) and (2) of BCG09 
for the definition of a and e), which does not affect the ac- 
creting object's structure and evolution, the intrinsic lumi- 
nosity of the protostar or brown dwarf is given by 




As mentioned above, we assume in the present calculations 
that the accreted matter reaches the stellar surface with a 
lower specific entropy than the object's internal one, which 
implies a — > in Eq. (0). As shown in BCG09, calcu- 
lations based on this assumption provide a consistent ex- 
planation of the luminosity spread in young cluster HRDs, 
while high values of a would predict overluminous objects 
during a short period of time (because of the short cor- 
responding Kelvin-Helmholtz timescale, see discussion in 
BCG09). Addition of fresh infalling deuterium and lithium 
is accounted for during the accretion process, with interstel- 
lar mass fractions [D]q = 2 x 10~ 5 and [Zi]o = 10 -9 , re- 
spectively. Depending on the accretion rate, these elements 
can start to be consumed at early stages in the proto-object 
interior during the burst accretion phases. The implemen- 
tation of the accretion process in the stellar evolution code 
has been validated (Gallardo 2007; Gallardo et al. 2009) by 
comparing our numerical results to those derived from the 
polytropic approach of Stahler (1988) and Hartmann et al. 
(1997) in the limit of validity of the polytropic model, i.e. 
for initial masses > 0.1 M©. 

2.2. Effect of episodic accretion on the internal structure 

As in BCG09, our accretion histories are based on the cal- 
culations of Vorobyov & Basu (2005), obtained from grav- 
itational instability in the accreting centrifugally supported 



disk. Starting from arbitrary initial masses Mj n ; t , we vary 
the number of bursts M, U rst, with accretion rate Mbm-st and 
duration Afburst, separated by quiescent phases of duration 
Af qu i et and accretion rate M qu j et . We explored a wide range 
of parameters, inferred from the results of Vorobyov & 
Basu (2005), with the aim of obtaining final LMSs within 
the mass range 0.1-1 M , in order to focus on the typical 
population used to derive the ages of star-forming regions 
and young clusters. Brown dwarf cases have been shown 
in BCG09. Our initial masses M- m \ x range from 1 Mj up to 
0.2 M Q . The typical burst number, accretion rates, and du- 
ration lie in the ranges A^m-st = 2 - 30, Mbm-st = 10 -4 - 5 
10~ 4 M yr _1 , and Af bu rst = 100 - 500 yr, respectively, while 
quiescent phases last between Af qu ; et = 10 3 and 10 4 yr with 
accretion rates M < 10~ 6 M o yr _1 . As mentioned in BCG09, 
accretion rates during the quiescent phase below this value 
have no significant impact on the final internal structure of 
the accreting object. For the sake of simplicity, we have thus 
adopted M qu ; et =0 in the present analysis, however, this ac- 
cretion luminosity must be taken into account when compar- 
ing the final luminosity with observations of proto-objects 
during the embedded phase. The next two sections illus- 
trate the typical results for two examples at each end of the 
presently studied mass range. 

2.2.1 . Sequences leading to the formation of a 0.1 M 
star 

Various evolutionary sequences of episodic accretion, with 
different initial masses, leading to a final 0. 1 M LMS are 
portrayed in Fig. Q] As mentioned in BCG09, high accre- 
tion rates yield significantly more compact structures, i.e., 
smaller radii compared to those of non accreting objects of 
same mass and age. This contraction comes from the in- 
crease in gravitational energy as mass is added, yielding 
higher central pressures and temperatures compared with 
the non accreting case. If the accretion timescale, = 
M t /M, is higher than the Kelvin-Helmholtz timescale, 
tkh = GMl/(R t L*), the object can relax to the radius 
it would have in the absence of accretion. In contrast, if 
t^j < tkh, which is the case for M burst = 10~ 4 - 5 10~ 4 
M B yr , the structure has no time to adjust to the incoming 
mass and energy, and the radius remains smaller than the 
non accreting counterpart of same mass and age, as illus- 
trated in Fig. Q] The more compact structure of the accret- 
ing object yields higher central temperature compared to the 
non accreting counterpart. Figure Q] highlights the main ef- 
fects of episodic accretion on proto low-mass objects, which 
can be summarized as follows. 

- For a given initial mass M m j t , the higher Mbm-st dur- 
ing the burst phases, the smaller the radius of the ob- 
ject at the end of the accretion phase, hence the higher 
its central temperature. This can be seen by comparing 
the long-dashed blue curve (Mj n ; t = 1 Mj up , Mburst = 
5 lO~ 4 M yr~') and the dash-dotted magenta curve (Mj n ; t 
= 1 M Jup , M burst = lO^yr 1 ). 

- For a given burst accretion rate Mbm-st, the lower the 
initial mass M m \ t , the greater the impact on the inter- 
nal structure, i.e, the smaller the radius and the higher 
the central temperature at the end of the accretion his- 
tory. This is illustrated by the comparison between the 
dash-dotted magenta curve (Mj n i t = 1 Mj up , Mb urs t 

1 "' If yr ' ) with the dotted red curve (M init = 10 M Jup , 
M burst = lO^Moyr 1 ). 



Baraffe & Chabrier: Effect of episodic accretion on internal structure and lithium depletion 



3 



It is worth stressing that the sequences starting with 
Mj n i t =l Mj up in Fig.[T]reach high enough central tempera- 
tures during the burst accretion episodes to completely burn 
deuterium, including the extra accreted fresh one. This ex- 
tra supply of nuclear energy is not sufficient to balance the 
increase in gravitational energy, hence the contraction of 
the object induced by the added matter. Deuterium burning 
starts when the central temperature T c reaches ~ 10 6 K. For 
the two sequences starting with Mj n i t =l Mj up and Mb urs t = 
10~ 4 M G yr _I and Mburst = 5 lO~ 4 M yr~', this corresponds 
to central densities of p c = 20 and 30 g ctrT 3 , respectively. 
This stage is reached when both sequences reach a mass 
of ~ O.O3M . For the sequence starting with M ln j t =10 Mj up 
and Mbmst = 10 _4 M o yr _I , the same central temperature is 
reached near the end of the accretion phase, for a central 
density of ~ 2 g ctrT 3 and a total mass ~ O.O8M . For 
this sequence, a significant fraction of deuterium (~ 65% 
of the initial mass fraction in the fully convective interior) 
is still present at the end of the accretion phase, generating 
an important source of nuclear energy which produces the 
increase in R and then decrease in T c between 10 4 yr and 
10 5 yr, as shown in Fig.Q] 

Finally, Fig.Q]clearly illustrates, for the sequences start- 
ing with 1 Mj up , the strong differences in radius and cen- 
tral temperature after a few Myr, up to 30 Myr, between a 
0.1 M star produced by episodic accretion with Mb urs t > 
lO~ 4 M yr~' and the non accreting counterpart. We also 
stress that, for initial masses M; n j t > O.O5M and burst accre- 
tion rates Mb urs t < 5 x lO~ 4 M yr~', accretion history hardly 
affects the internal structure and location in the HRD after 
~ 1 Myr, producing objects of 0. 1 M with properties simi- 
lar to those predicted by standard (non accreting) pre-main 
sequence models at ages > 1 Myr. 

2.2.2. Sequences leading to the formation of a 1 M Q 
star 

Figure [2] displays the same analysis for episodic accretion 
sequences producing 1 M objects. The figure portrays the 
results with Mburst = 5 x 10~ 4 M o yr _1 , which yields the 
strongest and most interesting effects on the structure and on 
lithium depletion (see §0. As found previously, the lower 
the initial mass, the greater the effect on the structure for 
a given accretion rate. For initial masses Mj n i t < 0.1 M 
and such high burst accretion rates, the figure illustrates the 
strong departure in radius and central temperature of these 
models from the non accreting sequence up to 20-30 Myr. 
For higher initial masses and Mb urs t < 5 x 10~ 4 M o yr _1 , the 
accretion history is found not to significantly alter the prop- 
erties of the newly formed 1 M star, producing an object 
with a position in the HRD and internal structure similar to 
the predictions of a non accreting evolutionary sequence at 
ages > 1 Myr. 

3. Effect on the size of the convective envelope 
and on lithium depletion 

3.1. The case of fully convective stars 

The significantly higher central temperatures induced by 
episodic accretion affect the timescale of lithium deple- 
tion, leading to faster lithium depletion compared with a 
non accreting object of the same age, as illustrated in Fig. 
[3] for sequences producing 0.1 M stars. Initial M m \ t = 
lMj U p protostars experiencing burst accretion rates Mb urs t = 




log t 



Fig. 1. Evolution of the radius (upper panel) and of the cen- 
tral temperature (lower panel) as a function of time (in yr) 
for models with episodic accretion and reaching a final mass 
of 0.1 M (see ^2T27Tb . Long-dash (blue): M init = 1 M Jup , 
Mburst = 5 10" 4 M Q yr 1 , AW =2; dash-dot (magenta): M init 
= 1 M Jup , Mb urs t = 10 -4 M Q yr-\ AW =10; dot (red): M M 
= 10 Mj up , Mburst = lO~ 4 M yr _1 , AWt =9. All calculations 
are done with Afb urs t=100 yr and Af qu i e t=1000 yr. The solid 
line (black) corresponds to the evolution of a non accreting 
0. 1 M star. 



5 10~ 4 M G yr _1 produce a 0. 1 M object that entirely depletes 
its lithium content within about 10 Myr, whereas for the 
non accreting counterpart, complete lithium depletion takes 
more than 50 Myr. As clearly illustrated in Fig. [3] depend- 
ing on the initial mass and the burst accretion rate, different 
episodic acccretion histories can produce objects with the 
same mass, say 0. 1 M , at a same age, say < 50 Myr, ex- 
hibiting different levels of lithium depletion. Similar effects 
are found for episodic accretion sequences producing final 
objects within the entire characteristic domain of fully con- 
vective stars, M < O.35M . 



3.2. The case of stars that develop a radiative core 

If the central temperature exceeds ~ 2-3 10 6 K, a radiative 
core develops, because of the opacitiy decrease after the last 
"/<--bump" due to metals (C, O, Ne and Fe, see Rogers & 
Iglesias 1992), as explained in Chabrier & Baraffe (1997, 
see their §3.2 and their Fig. 9). The exact temperature at 
which this occurs depends on the density, given the sensitiv- 
ity of the opacities to density in this temperature range (see 
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Fig. 2. Same as Fig.Q]for models reaching a final mass 1 M 
(see $FTR . Long-dash (blue): M init = 10 M Jup , AW =20; 
dash-dot (magenta): M- m \ t = 30 Mj up , AW-st =20; dot (red): 
Mnit = 0.1 M Q , A^burst =18. All calculations are done with 
Mburst = 5 lO^Moyr 1 , Af burst =100 yr and Af quiet =1000 yr. 
The solid line (black) corresponds to the evolution of a non 
accreting 1 M & star. 
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Fig. 3. Evolution of the lithium abundance (divided by the 
initial Li abundance) as a function of time for the same mod- 
els producing a 0.1 M low-mass star as in Fig.Q] 



depletion, yielding very different final lithium abundances 
in the convective envelope, as illustrated in Fig. [4] 

The main results of this section can be summarized as 
follows. Depending on the initial seed mass and the burst 
rates, different episodic accretion histories can produce 1 
M Q star models with different sizes of the convective enve- 
lope at the same age, for ages < 30 Myr. The models can 
have different surface Li abundances even after the models 
have converged toward the same structure, i.e. after ~ 30 
Myr in the particular case portrayed in Fig. [4] Similar ef- 
fects are found for models with episodic accretion produc- 
ing partly convective objects in the mass range 0.35 -1 M . 



Fig. 2a of Rogers & Iglesias 1992). The higher the density, 
the higher the temperature required for the radiative core 
to develop. Since this range of temperatures is also char- 
acteristic of the temperature required for Li nuclear fusion, 
the central temperature and density at which the radiative 
core develops will determine the temperature at the bottom 
of the convective envelope and thus the level of Li deple- 
tion in the convective envelope. For sequences producing 1 
M stars, Fig. [4] shows that accretion history has a strong 
impact on (i) the age for the onset of a radiative core, (ii) 
the mass of the convective envelope at ages < 30 Myr, and 
(iii) the Li abundance in the convective envelope. The higher 
central temperature reached by sequences with strong burst 
accretion rates (see Fig. [2]) results in a radiative core that de- 
velops earlier in time. For the sequence starting from M; n j t 
= 10 Mj up , its more compact structure yields significantly 
higher temperatures at the bottom of the convective enve- 
lope, with a maximum of ~ 7 10 6 K, resulting in complete 
Li depletion in the convective envelope at ages < 1 Myr. In 
comparison, for the accreting sequence starting with M- m \ t 
= 0.1 M , the maximum temperature reached at the bottom 
of the convective envelope is ~ 3.5 10 6 K, which is compa- 
rable to the value found in the non accreting sequence (see 
Fig. 3). These differences in temperature are crucial for Li 



4. Discussion 

4. 1 . Summary of the results 

We have shown in the present work that an early proto- 
star/BD accretion history based on episodes of short, in- 
tense bursts of accretion with typical accretion rates Mb urs t = 
10~ 4 - 5 lO _4 M yr _1 , as obtained in 2D hydrodynami- 
cal simulations of gravitationally unstable accretion disks 
(Vorobyov & Basu 2005), can affect the internal structure of 
low-mass objects even after several Myr, up to a few times 
10 Myr. The main results obtained in the present analysis 
can be summarized as follows. 

- (i) Episodic accretion produces objects with smaller ra- 
dius and higher central temperature compared to the non 
accreting counterpart with the same mass at the same 
age. 

- (ii) Higher central temperatures at a given age and mass 
can significantly enhance lithium depletion in fully con- 
vective low-mass stars. 

- (iii) As a consequence of the hotter structure, a radia- 
tive core develops earlier in accreting objects with final 
masses M >O.35M . Episodic accretion can thus lead to 
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Fig. 4. Evolution of the mass of the radiative core divided 
by the total mass (upper panel) and of the surface lithium 
abundance divided by the initial Li abundance (lower panel) 
as a function of time for models producing a 1 M star. The 
curves refer to the same models as in Fig. [2] Note that Li is 
entirely depleted in the model with Mi n ; t =10 Mj up (long- 
dashed blue line) in less than 1 Myr. 



final objects with smaller convective envelopes than pre- 
dicted by standard non accreting models, at ages up to a 
few times 10 Myr. 
- (iv) The more compact and hotter structure of accret- 
ing models with M>O.35M increases the maximum 
temperature reached at the bottom of the convective en- 
velope, compared to the non accreting counterpart, in- 
creasing the level of lithium depletion in the convective 
envelope. 

We stress that the magnitude of all these effects strongly de- 
pends on the initial mass, the acretion burst rate, and the 
total number of bursts, but barely depends on the phases 
of quiescent accretion, as long as M qu i et <k 10~ 6 M o yr _1 . 
Within the range of parameters mentioned in 92.21 episodic 
accretion can produce final low-mass star models that de- 
part at various levels from standard non accreting models in 
terms of radius, luminosity (and thus position in the HRD, 
see BCG09), lithium abundance, and convective envelope 
mass, at a given similar age. Conversely, final objects is- 
sued from episodic accretion display observational signa- 
tures similar to non-accreting objects of same mass, but at a 
younger age, as shown in BCG09, the age difference again 
depending on the various accretion histories. 



4.2. Comparison to observations: anomalous lithium 
depletion in young clusters 

Several observational studies report discrepancies between 
ages inferred from evolutionary isochrones, on one hand, 
and from lithium depletion, on the other hand, for M to 
G type stars of various ages (e.g Yee & Jensen 2010). 
Observations also show lithium scatter in clusters and as- 
sociations (e.g da Silva et al. 2009; King et al. 2010) and 
young cluster members exhibiting severe, unexpected Li de- 
pletion for the inferred age of the cluster (e.g Kenyon et al. 
2005; Sacco et al. 2007). The present calculations provide 
a consistent solution to these puzzles, in particular provid- 
ing an explanation for the anomalous lithium depletion ob- 
served in some members of clusters that are a few Myr old 
(Kenyon et al. 2005; Sacco et al. 2007; Sacco et al. 2008; da 
Silva et al. 2009 ). The explanation suggested for these "in- 
terlopers", found for instance in cr Ori, A Ori, and the ONC, 
is to invoke a significant age spread within these clusters. 
However, at least in some cases (see below), the required 
age spread is significant, up to 10-20 Myr, a rather unlikely 
possibility in such young clusters. The present calculations, 
based on early evolution sequences involving accretion out- 
bursts, provide a more plausible explanation, as illustrated 
in Fig. [5] for the three weird objects discovered by Sacco 
et al. (2007) in cr Orionis. These objects show severe Li 
depletion, inconsistent with the age of the cluster (about 
5 Myr), but have radial velocities consistent with cluster 
membership. They are displayed in the upper panel of Fig. 
[5] Evolutionary sequences undergoing burst accretion with 
rates Mburst ~ 5 10 _4 M Q yr _1 and reaching final masses 0.3 
M Q and 0.7 M Q , respectively, reach the observed locations 
in the HRD at an age of ~ 5 Myr, consistent with the age of 
cr-Ori, whereas the same locations correspond to 10-15 Myr 
for the non accreting sequences of same mass, as shown in 
the figure. For the accreting sequence producing a 0.7 M 
object, the middle panel of Fig.[5]shows that a radiative core 
develops much earlier compared to the corresponding non 
accreting sequence. The lower panel illustrates the much 
faster Li depletion for the accreting sequences, as discussed 
in $3] with complete destruction occurring within less than 
1 Myr for the sequence producing a 0.7 M object. The 
possibility for the accretion burst models to explain consis- 
tently both the position in the HRD and the observed severe 
lithium depletion for these objects at the cluster age without 
particular fine tuning of the accretion parameters provides 
a strong support for this scenario. Although the high level 
of Li depletion in the three objects of Sacco et al. (2007) 
has recently been questioned by Caballero (2010), similar 
anomalous Li-depletions have been reported in other clus- 
ters by other groups (see references above), which can also 
be explained by episodic accretion. Our results thus rein- 
force the word of caution of da Silva et al. (2010) regarding 
the elimination of stars as cluster members based only on Li 
abundances. The present analysis demonstrates that lithium 
is not a reliable age indicator, as its fate strongly depends on 
the past accretion history experienced by the star! 

The same episodic accretion calculations can explain the 
unexpected degree of lithium depletion of stars with mass > 
0.5 M in the young cluster IC 4665 (Jeffries et al. 2009), or 
the spread in Li abundance in cool dwarfs of about a solar 
mass reported in the Pleiades (King et al. 2010), although 
other scenarios, based either on the star rotational history 
(Bouvier 2008) or stellar activity (King et al. 2010), have 
been suggested. 
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Fig. 5. Comparison of accreting models with observations of Sacco et al. (2007). Upper panel: HR diagram. The two long- 
dashed (black) curves are the 1 Myr and 10 Myr isochrones of Baraffe et al. (1998) for non accreting models. The black 
points are the objects of Sacco et al. (2007). Middle panel: radiative core mass versus time for models producing a 0.7 
M Q star. Lower panel: surface Li abundance versus time. The solid red curve is indistinguishable from the X-axis (i.e. 
[Li] = 0): accreting sequence producing an object with final mass Mf=0.7 M with Mi n ; t = 20 Mj up , Nbmst = 14; Dash- 
dotted blue: accreting sequence producing an object with final mass Mf=0.3 M with M- m \ x = 10 Mj up , AW- s t =6. In both 
cases, Mburst = 5 10 _4 M G yr _1 , Afb urs t=100 yr and Af qu j et =1000 yr; Dashed red: non accreting 0.7 M model; Dotted blue: 
non accreting 0.3 M model. 



4.3. Abundances in planet host stars 

Another interesting implication of our calculations and the 
inferred link between episodic accretion, lithium depletion, 
and the size of the convective envelope in young stars is 
the connection with the presence of close-in giant planets. 
Zhu et al. (2010b) suggest that the presence of mass clumps 
in some specific regions of circumstellar disks, which may 
lead eventually to non steady accretion when the disk be- 



comes unstable (Vorobyov & Basu 2005), could play an im- 
portant role in the process of planet formation. Following up 
along these lines, we can speculate that the presence of giant 
planets, which require massive protoplanetary disks, more 
prone to gravitational instability, could be linked to the in- 
tensity of episodic accretion outburts. If this suggestion hap- 
pens to be true, the present calculations provide an explana- 
tion for the claimed higher lithium depletion in planet host 
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stars (Israelian et al. 2009; Sousa et al. 2010): since these 
planet host stars will have undergone more intense accre- 
tion bursts, they will experience more intense Li depletion 
than similar planetless stars whose presumably less massive 
disks are less prone to violent accretion-burst episodes. 

Furthermore, as shown in this paper, episodic accre- 
tion can produce objects whith smaller convective envelopes 
than the ones obtained with standard non accreting mod- 
els at ages < 10 Myr (see Figs.|4]and|5]), the characteristic 
maximum lifetime of protoplanetary disks. Consequently, 
we can expect the signature of (both terrestrial and gaseous) 
planet formation to leave a greater imprint on stars that 
have undergone phases of intense accretion bursts. Although 
speculative, this suggestion provides a consistent explana- 
tion for the recently claimed peculiar chemical abundance 
of the Sun, with about 20% depletion of refractory ele- 
ments relative to volatile elements, compared with solar- 
type analogs harboring close-in giant planets (Melendez et 
al. 2009; Ramirez et al. 2009). According to these authors, 
this chemical peculiarity seems to be common among Sun- 
like stars with no detected giant planets and is said to reflect 
the signature of the formation of terrestrial planets, which 
trap the refractories, whereas the dust-depleted gas is ac- 
creted by the central star (Ramirez et al. 2009; Melendez et 
al. 2009). For such a suggestion to be viable, and the sig- 
nature of terrestrial planet formation to be imprinted in the 
Sun photospheric composition, the Sun's convective enve- 
lope, of mass M env , must have been much smaller during 
the lifetime of the accretion disk, i.e. during the first ~ 10 
Myr, than predicted by standard pre-main sequence models. 
Indeed, let X'. and X, be the mass fractions of element i in 
the accreting material and in the originally unpolluted star 
convective envelope, respectively. The new abundance X'.' 
of such element in the convective envelope becomes 



x': 



-^-d+ A 

l+a X, J 



(5) 



where a = M acc /M env is the ratio of the accreted mass to the 
star convective envelope mass. The ratio of two elements, 
say i =refractories and j =volatiles, in the convective enve- 
lope after accretion of disk material becomes 



X" 

J 



1 x 'i 
Xj i + "t 

v. X' ' 



(6) 



In the presence of terrestrial planets, we expect X\ < Xj 
and Xj > Xj, since most refractories remain trapped in the 
planet, and thus 

(7) 



X" 



X, 



The larger a, hence the smaller the convective envelope, the 
larger the inequality Eq.©, hence the larger the indirect sig- 
nature of planet formation in the star's convective envelope. 

To obtain such a result, the aforementioned authors re- 
fer to the results of Wuchterl & Klessen (2001), who ex- 
plored the formation and early evolution of a 1 M star, and 
found that the early Sun was never fully convective. These 
calculations, however, are incorrect. First of all, as dis- 
cussed in Baraffe et al. (2003), the calculations of Wuchterl 
& Klessen (2001) yield an evolutionary track for a 1 M 
star that is much too hot, by ~ 700 - 800 K in T e g at 
a given L, compared to observed young binary systems 
with (dynamically determined) masses around a solar mass. 



Second of all, the unrealistic assumption of spherical accre- 
tion in the Wuchterl & Klessen calculations leads to signif- 
icantly hotter inner structures, which will favor the growth 
of a radiative core, compared with more realistic 3D rota- 
tional collapse calculations (Chabrier et al. 2007, Fig. 3). 
This assumption is thus certainly responsible for the above- 
mentioned large overestimate of the effective temperature at 
young ages. On the other hand, the effect of episodic ac- 
cretion on the mass of the convective envelope illustrated in 
Fig. |4] provides a more plausible explanation for this prob- 
lem, and suggests that the early Sun could have undergone 
a phase of strong accreting bursts. 

Although interesting, however, this interpretation of the 
peculiar refractory to volatile abundance ratio in the Sun and 
other solar twins without detected close-in giant planets as 
a consequence of terrestrial planet formation is challenged 
by the recent analysis of Chavero et al. (2010), based on 
four CoRoT planet host stars. This work suggests that alter- 
ation of abundances of refractories with respect to volatiles 
in stellar atmospheres may simply result from condensation 
processes in the accretion disk and from accretion of such 
altered gas onto the star. Independent of the source of this 
peculiar abundance pattern, the signature of such conden- 
sation/accretion processes will still leave a greater imprint 
if the star convective envelope is considerably smaller than 
usually expected for solar type stars at ages <10 Myr. Also, 
the analysis of Chavero et al. (2010) reveals no obvious 
correlation between abundances and condensation temper- 
atures in the CoRoT stars, indicating no sign of overabun- 
dance of volatiles relative to refractories in stars harboring 
close-in giant planets compared to the Sun, contradicting the 
observations of Melendez et al. (2009) and Ramirez et al. 
(2009). The peculiar abundance ratio determinations in the 
Sun and solar analogs without detected giant planets must 
thus be confirmed by further studies. 

In conclusion, anomalous abundance ratios of refracto- 
ries to volatiles in stars may not necessarily reflect the signa- 
ture of terrestrial planet formation but may also result from 
condensation processes in the disk. A detectable signature 
of these abundance anomalies, however, requires a smaller 
convective envelope during typical disk lifetimes than pre- 
dicted by standard models. Strong accretion outburts do pro- 
duce these favorable conditions. 



5. Conclusion 

We have shown here that non steady accretion can strongly 
affect the internal mechanical (radius) and thermal (temper- 
ature) structures, hence the observational signatures (L, T e s) 
of low-mass stars (and brown dwarfs), possibly up to ages of 
a few tens of Myr. The magnitude of the effects strongly de- 
pends on the forming object's initial core mass, the number 
of bursts, and the intensity of the burst accretion episodes. 
As already highlighted in BCG09, episodic accretion can 
produce models of a given final mass that to various ex- 
tents depart from standard non accreting models in terms of 
position in the HRD, Li abundance, and convective enve- 
lope mass. Since the mechanisms responsible for the accret- 
ing disk formation and possible outbursts depend on many 
environmental parameters (e.g. mass, temperature, angular 
momentum, magnetic field of the parent collapsing cloud, 
and thermal properties of the disk), it is certainly possible 
to produce various populations of stars at a similar age in 
young clusters and associations that depart at various lev- 
els from predictions of standard pre-main sequence mod- 
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els. Therefore, our scenario does not systematically predict 
a broad spread among these stellar properties/signatures, be- 
cause the level of departure entirely depends on the star or 
cluster properties and accretion histories. Further observa- 
tional effort is needed to explore these issues in more de- 
tail and to try to have better constraints on the processes of 
non steady accretion during the early stages of evolution. 
Among the multiple consequences of early episodic accre- 
tion, we emphasize the limited (to say the least) reliability of 
Li depletion as an age or cluster membership indicator. We 
urge observers to re-analyze young cluster populations on 
the basis of the present results, since several genuine cluster 
members may have been eliminated because of their unex- 
pectedly low Li abundance. Any information about the frac- 
tion of such Li depleted objects, with membership unam- 
biguously assessed by radial velocity, could provide crucial 
insight into the occurrence and the intensity of non steady 
accretion. 

We also suggest that intense bursts of episodic accretion 
may have some link to planet formation, ;/ planet forma- 
tion preferentially occurs in massive disks, which are more 
prone to instability. We speculate that stars harboring plan- 
ets, including our Sun, may have experienced burst accre- 
tion episodes during their youth, which then led to a hot- 
ter thermal structure, thus to faster lithium depletion and a 
smaller outer convective zone than conventionally admitted. 
This provides a natural explanation for the higher lithium 
depletion observed in stars harboring at least giant planets. 
This also provides a plausible, although admittedly specu- 
lative, explanation for the claimed underabundance of re- 
fractory elements in the atmosphere of stars with no de- 
tected close-in giant planets, as a result of the formation 
of terrestrial planets, which trap large amounts of refrac- 
tory elements. In any event, the multiple consequences of 
episodic accretion that we characterize in the present work 
should stimulate deeper theoretical and observational in- 
vestigations, in order to better understand the twilight zone 
characteristic of (proto)star early evolution. 
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